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Abstract. The combination of ground-based astrometric com- 
pilation catalogues, such as the FK5 or the GC, with the re- 
sults of the ESA Astrometric Satellite HIPPARCOS produces 
for many thousands of stars proper motions which are signifi- 
cantly more accurate than the proper motions derived from the 
HIPPARCOS observations alone. In the combination of the ba- 
sic FK5 with the HIPPARCOS Catalogue (i.e., in the FK6), the 
gain in accuracy is about a factor of two for the proper mo- 
tions of single stars. The use of the GC still improves the accu- 
racy of the proper motions by a factor of about 1 .2 . We derive 
and describe in detail how to combine a ground-based com- 
pilation catalogue with HIPPARCOS. Our analytic approach 
is helpful for understanding the principles of the combina- 
tion method. In real applications we use a numerical approach 
which avoids some (minor) approximations made in the ana- 
lytic approach. We give a numerical example of our combina- 
tion method and present an overall error budget for the combi- 
nation of the ground-based data for the basic FK5 stars and for 
the GC stars with the HIPPARCOS observations. 

In the present paper we describe the 'single-star mode' of 
our combination method. This mode is appropriate for truly 
single stars or for stars which can be treated like single stars. 
The specific handling of binaries will be discussed in subse- 
quent papers. 
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1. Introduction 

The HIPPARCOS astrometric satellite has provided very ac- 
curate positions, proper motions, and parallaxes for more than 
118 000 stars (ESA 1997). Does this mean that the long se- 
ries of ground-based astrometric observations of all these stars 
are now merely of historical value ? This is certainly not the 
case: It can be shown that the combination of the HIPPAR- 
COS data with ground-based results is providing for many 
thousands of stars individual proper motions which are signif- 
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icantly more accurate than the HIPPARCOS proper motions 
themselves (Wielen 1988, Wielen et al. 1998, 1999). 

There are two main reasons why ground-based observa- 
tions can improve the HIPPARCOS results: (1) For many of 
the brighter stars, the ground-based observations cover a period 
of time of more than two centuries. The measuring accuracy 
of these ground-based positions is such that they allow, espe- 
cially in combination with the HIPPARCOS positions, to derive 
proper motions which have significantly smaller measuring er- 
rors than the HIPPARCOS proper motions. It is well-known 
that ground-based observations suffer from considerable sys- 
tematic errors. However, it is just the HIPPARCOS Catalogue 
itself which allows us to remove these systematic errors from 
the ground-based data to such a high degree that we can use 
with much confidence the corrected ground-based results for 
the direct combination with HIPPARCOS. (2) The HIPPAR- 
COS proper motions are derived from observations which have 
been carried out within about three years only. Since many 
of the HIPPARCOS stars are undetected astrometric binaries, 
these nearly 'instantaneously' measured HIPPARCOS proper 
motions can differ significantly from 'time-averaged', 'mean' 
proper motions. We have called the difference between the in- 
stantaneous HIPPARCOS proper motion and the mean motion 
of a star the 'cosmic error' of the HIPPARCOS proper motion 
(Wielen 1997). In many cases these cosmic errors are signif- 
icantly larger than the measuring errors of the instantaneous 
HIPPARCOS proper motions (Wielen 1995a, b, Wielen et al. 
1997, 1998, 1999). The ground-based observations, with their 
long observational history, are often already providing 'mean' 
proper motions. The ground-based data allow us therefore to 
identify and to correct partially the cosmic errors in the HIP- 
PARCOS proper motions. 

Having argued that a combination of ground-based obser- 
vations with the HIPPARCOS data is useful for many stars, we 
now provide a method for carrying out such a combination. The 
best way would be to combine the individual observational cat- 
alogues of ground-based positions directly with HIPPARCOS. 
However, there exists more than 2000 of such observational 
catalogues. While we have collected most of these catalogues 
in our astrometric data base ARIGFH (described in Wielen 
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(1998)) and hope to use them individually in the future, we 
have to rely for the moment on a few compilation catalogues. 
Such a compilation catalogue is usually based on a large num- 
ber of observational catalogues and provides as a 'summary' 
of these individual catalogues a derived proper motion and a 
mean position at a central epoch. 

The most accurate of the ground-based compilation cata- 
logues is the first part of the FK5 (Fricke et al. 1988). It contains 
the 1535 basic fundamental stars. The second part of the FK5 
(Fricke et al. 1991) lists 3117 additional stars in the bright and 
faint extension. A compilation of the remaining FK4Sup stars 
(Schwan et al. 1993) provides data for further 995 stars with a 
good observational history. The combination of the basic FK5 
stars with HIPPARCOS data produces the FK6. The proper mo- 
tions of the FK6 are of unprecedented accuracy (Wielen et al. 
1998,1999, and Sect. 6). The accuracy of these FK6 proper mo- 
tions rests to a high degree on the older observations. This is not 
only to be expected on theoretical grounds, but can be proved 
empirically also (Wielen et al. 1998) by combining HIPPAR- 
COS with an older fundamental catalogue such as the FK3 
(Kopff 1937, 1938). Because the old observations carry most of 
the weight in a combination with HIPPARCOS, the GC (Boss 
et al. 1937) is also relevant for our purpose. The GC is a careful 
compilation of ground-based observations from the 18th cen- 
tury until about 1930. The main advantage of the GC is its 
large number of stars. About 29 700 of the 33 342 GC stars 
have been observed by HIPPARCOS and can therefore be used 
in our method. The typical accuracy obtained for a GC star is, 
of course, lower than for an FK6 star (Wielen et al. 1998, 1999). 

Instead of using a compilation catalogue, we may rely on 
one old observational catalogue only, such as the Astrographic 
Catalogue (AC). However, the formal accuracy of the proper 
motions obtained by a combination of the AC with the HIP- 
PARCOS data is in most cases only marginally better than that 
of the HIPPARCOS proper motions alone. On the other hand, 
a comparison of a proper motion derived from an AC position 
and a HIPPARCOS position with the quasi-instantaneous HIP- 
PARCOS proper motion [in can reveal the occurence of a large 
value of the cosmic error in [in for the star under consideration. 

Before we can combine the ground-based data with HIP- 
PARCOS, the systematic errors of the ground-based positions 
and proper motions have to be determined and removed. For the 
determination of these systematic errors we use methods devel- 
oped for the construction of the FK5 (Bien et al. 1978). These 
methods do not only provide the systematic errors themselves 
but also the local uncertainty of the correction, i.e. the mean 
'measuring' error of the systematic correction for a given star. 
As the reference system, we use the HIPPARCOS system. The 
combined catalogue is therefore always on the HIPPARCOS 
system, i.e. on the ICRS. For our purpose, we can consider the 
HIPPARCOS data as free from systematic errors. A possible 
slight rotation of the HIPPARCOS system with respect to an 
inertial, extragalactic system does not affect our method (but 
of course our results). 

In the following sections, we assume that the angular coor- 
dinates of a star, a* = a cos 8 and 5, change linearly in time 



t. In practice, the non-linear motion of a star in a* and 8 be- 
cause of spherical effects and the foreshortening effect has to be 
taken into account. We do this already in the determination of 
the systematic errors between the ground-based data and HIP- 
PARCOS. Later we do not use the full values of ground-based 
positions and proper motions, but only their differential values 
with respect to the HIPPARCOS solution at the proper epoch. 
If we use these differential values only, the linear approxima- 
tion is then fully sufficient for our purpose. 

We shall now present the mathematical details of our com- 
bination method. We offer two approaches: an approximate an- 
alytic one and a rigorous numerical one. The analytic approach 
has the advantage to show more clearly the basic principles. It 
neglects, however, some of the correlations which occur among 
the HIPPARCOS results. In practical applications, we use the 
numerical approach which takes all the correlation among the 
HIPPARCOS results into account. In all cases, we assume that 
the HIPPARCOS data are completely uncorrelated with the 
ground-based data. This is certainly true for all of our stars, 
although the HIPPARCOS Input Catalogue, used as a first ap- 
proximation in the HIPPARCOS data reduction procedure, had 
to rely on the ground-based observations. 

2. Analytic approach 

We assume that two astrometric catalogues are available. They 
are identified by the indices 1 and 2. Examples are the FK5 and 
HIPPARCOS. For the combined catalogue, e.g. the FK6, we 
use the index C. Each of the two basic catalogues (i = 1, 2) pro- 
vides for all the stars under consideration a position X{ (Tj) and 
a proper motion [n at a central epoch Tj for two uncorrelated 
coordinate components, e.g. for a* and 8. In addition, both cat- 
alogues provide mean errors of Xi(Ti) and in which we denote 
by e X)j and e^. 

In the case of a catalogue which has been reduced to the 
HIPPARCOS system, the errors s Xt i and e^.i have to take into 
account also the local uncertainty of the systematic corrections, 



e.g. by using e% l tot 



+ et 



where the 'individ- 



ual' error z x ,i,ind is the random measuring error given in the 
catalogue. We should emphasize that e x ^ tSys is not the system- 
atic correction itself but only its uncertainty. 

The central epoch Tj is chosen such that the correlation be- 
tween Xi(Ti) and [n is zero at t = T{. Most of the ground- 
based catalogues give the central epoch T, explicitely for each 
star, separately for a and 8. For HIPPARCOS (index 2) we have 
to derive T 2 = Tn,ind (separately for a* and 8) from the data 
given in the catalogue: 



T 2 =T, 



H.ind 



T. 



H 



Pxn,H (Th) Sx,H (Th) 



(1) 



£fi,H (Th) 

where Tjj — 1991.25 is the overall reference epoch of the HIP- 
PARCOS Catalogue, s Xt H (Th) the mean error of x H (Th), 
e^.h the mean error of ^h, and p x ^H (Th) the correlation 
coefficient between xh(Th) and \ih- The position x 2 (T 2 ) is 
derived from 



x 2 (T 2 ) = x H (T H ) + ^h(T 2 - T H ) ■ 



(2) 
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The mean error of X2 (T 2 ) is given by 



catalogue i. The normal equations of catalogue i are then 



+ 1p x ^ H {T H )e XtH {T H )e^ H {T 2 -T H ) . (3) 

The proper motion has not to be recalculated, i.e. ijl 2 = I^h and 
£fi,2 = £fi,H- The HIPPARCOS position and proper motion, 
2-2(^2) and 112, are now uncorrelated, i.e. p x ^,2 

(T 2 ) = 0. In 

the analytic approach, we assume that the values of £2(22) and 
[12 between a* and 8 are also uncorrelated. This is not strictly 
true because our choice of T 2 (separately for a* and 8) forces 
only one correlation coefficient to vanish. All the other correla- 
tion coefficients of the HIPPARCOS data remain finite (albeit 
small) in general. Practical experience shows that this approxi- 
mation is quite good for most stars (see Sect. 5). 

We are now prepared to carry out the combination of the 
two catalogues. The basic idea of the combination is to recon- 
struct the normal equations of the least-square solutions from 
which the two catalogues have been derived. This can be done 
by using only the data actually given in the two catalogues. 
These two sets of normal equations are then added together to 
provide the normal equations for the combined catalogue C. 
Our method is a generalisation of methods described by Kopff 
et al. (1964), Eichhorn (1974), and other authors. First partial 
results were presented by Wielen (1988). 

The method of least squares determines the 'best' solution 
x so i(T re f) and /i so i for the position of a star at a reference 
epoch T re f and for its proper motion from a time series of ob- 
served positions, 



x bs(t) = x so i(T ref ) + fi so i{t - T ref ) + v obs (t) , 



(4) 



by the condition that the weighted sum of the squared residuals 
v bs {t) should be a minimum: 



[P v lbs] = min . 



(5) 



We use here the classical convention that the brackets [...] 
imply a summation over the observations, p is the individual 
weight of x bs(t). (We hope that the reader is not confused by 
the fact that we use the same letter p for two different quan- 
tities in order not to change familiar notations: In this Sect. 2, 
we denote by p the weights of observations. In all the other sec- 
tions, p is the stellar parallax.) The condition (5) is fulfilled if 
we solve the normal equations for x so i(T re f) and p so i\ 

[p] x so i(T ref ) + \p(t - T ref )} ij, so i 

= \px bs(t)] , (6) 

\p(t - T re f)} X so l(T re f) + [p(t - T ref ) 2 ] (J, 3ol 

= \px obs (t)(t-T ref )] . (7) 

We now specialize the Eqs. (6) and (7) for each of the basic cat- 
alogues (i = 1, 2). For this purpose, we choose individual refer- 
ence epochs T re f for each catalogue (T re f = Tj) such that 



\p{t-Ti)]i = 0. 



The index i at the right bracket indicates that the sum should 
be extended over all the observations which have been used in 



[p(t - T t ) 2 ] Mj 



[px obs (t)] i 
[px obs (t)(t - 



(9) 
(10) 



The choice (8) of the central epoch Tj for T re f has the conse- 
quence that Xi(Ti) and in are not correlated, as indicated by 
the vanishing diagonal terms which leads to decoupled equa- 
tions for Xi(Ti) and m. 

The mean errors of Xi(Ti) and £ Mj i of in are given by: 



Mi ' 



(11) 



(12) 



where £0 is the error of unit weight. It will turn out that we do 
not have to know £0 or the adopted system of weights in each of 
the basic catalogues. We rely only on the fact that mean errors 
were correctly calculated for the final results given in each of 
the catalogues. 

The four quantities in the brackets of Eqs. (9) and (10) are 
now redetermined in the following. From Eqs. (11) and (12), 
we obtain 



(13) 



- X,l 

-2 



W-Ti)\ = d4) 

Since the catalogue values of ^(Tj) and Hi solve the Eqs. (9) 
and (10), we can reconstruct the right-hand sides of the normal 
equations with the help of Eqs. (13) and (14): 



[pX obs (t)] l 
IpXobsWit-Ti)^ 



2 Xi(Ti) , 



F 2 



(15) 
(16) 



The normal equations for the combined catalogue C have 
the same form as Eqs. (6) and (7). In order to simplify these 
normal equations, we replace T re f by the central epoch Tc of 
the combined catalogue C. Tc is given by the condition: 



[ P (t - T c )] c = 



(17) 



Tc can be determined from known quantities by using Eqs. (8): 

[ P (t - T C )] C = W - Tc)}, + [p(t - T C )) 2 
= [p(t-T 1 )} 1 + [p(T 1 -T c )} 1 
+ [p(t - T 2 )} 2 + [p(T 2 - T C )] 2 

= [p} 1 (T 1 -T c ) + [p} 2 (T 2 -Tc)=0 . (18) 
Solving for Tc and inserting Eqs. (13), we find 



(8) Tc - 



p- 2 +F~ 2 



Mi Ti + [ P ] 2 T 2 _ e x l T\ + e x 2 2 T 2 
w x ,i T\ + w x , 2 T 2 

W x .\ + W Xi 2 



(19) 
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The central epoch Tc is therefore the weighted mean of the 
epochs T\ and T 2 . As weights for T i7 we have to use the weight 



W x .i = 



(20) 



~x,l 



of the position Xi(Ti). Using Tc as the reference time, the nor- 
mal equations for the combined catalogue C are reduced to two 
independent equations for the unknowns xc{Tc) and \ic, the 
position and proper motion of the star in C: 



\p] c x c (T c ) = [px obs (t)] c , 
[p(t-T c ) 2 ] c [i C = [px bs{t){t-Tc)] c 



(21) 
(22) 



If we solve Eq. (21) for x c {T c ) and use Eqs. (9), (13), and 
(15), we find 



xc{Tc) 



[pX obs (t)] C _ \pX obs (t)] 1 + \pX obs (t)]r 



[p\ c Mi 

[p] x xi(Ti) + [p} 2 x 2 (T 2 ) 

Mi + [v\ 2 

w x .iXi{Tx) + w x , 2 x 2 (T 2 ) 

W x .l + W Xj2 



(23) 



with w xA given by Eq. (20). As for Tc, the position xc{Tc) is 
the weighted mean of X\{Ti) and x 2 {T 2 ) with w x .i as weights. 

In order to write hc also as a weighted average, it is help- 
ful to introduce a third proper motion, ho, as an auxiliary tool 
(Wielen 1988): 



Mo 



x 2 (T 2 ) - x^n) 



(24) 



The proper motion ho is based only on the central positions of 
the two basic catalogues (see also Fig. 1). Since £i(7i) and 
x 2 (T 2 ) are not correlated with the catalogue proper motions Hi 
and fi 2 , the new proper motion ho is also not correlated with 
Hi and [i 2 . The mean error e^o of ho is given by 



c x,l 



(T 2 - 7\)2 • 
The weights lo^j of the proper motions are 



w„ 



(25) 



(26) 



where the index i now runs from to 2. We now solve Eq. (22) 
for^ c : 



Mc = 



[px obs (t)(t - T c )] c 
W-Tc?] c 



(27) 



For the numerator of Eq.(27) we obtain by using Eqs. (15), 
(16), (19), (20), (24), (25), and (26): 

\px obs (t){t-T c )] c 

= [px obs (t){t - Tc)]! + [px obs (t){t - T c )} 2 
= [px obs (t){t - Ti^ + [px obs (t){t - T 2 )] 2 

+ (7\ - T c ) [px obs (t)} 1 + (T 2 - T c ) [px obs (t)} 2 



w x . 2 



+ ■ 



W x .l + W x>2 



{T 1 -T 2 )w x , 1 x 1 (T 1 ) 

(T2-T 1 )w x ,2X2(T 2 f) 



Wx.i + Wx,2 

el (w^ini +w^ 2 (J>2 

Wx.i w x ,2 (T 2 - Ti) 2 x 2 (T 2 ) -ari(Ti) 



+ ■ 



4 



W x ,l + Wx.2 
1 Ml + W^.2 M2 + — 



T 2 -Ti 
7\) 2 



) 



~'x,l 



+ e 



el (w M ,i Mi + w^.2 M2 + w M ,o Mo) • 



— Mo) 

x,2 ' 



(28) 



Similar to Eq. (28) we obtain for the denominator of Eq. (27) 
after some algebra 

[p(t Tc) 2 ] c = [p(t Tc) 2 ] 1 + [p(t Tc) 2 ] 2 
~ p {[t-T 1 ) + {T 1 -Tc)Y 



+ 



p((t-T 2 ) + (T 2 -T c )Y 



(29) 



(30) 



= [p(t-T 1 )\+[p(t-T 2 )*] 2 

+ (T 1 -Tc) 2 [p] 1 + (T 2 -T c ) 2 [p] 2 
= e 2 {w^+w^+Wpfl) ■ 

Our final result for He is therefore, using Eqs. (27) - (29), 

w^i Hi + uy 2 M2 + w^q /j, 
Mc = — — — • 

Wft,i + w^ 2 + w^ 

Equation (30) means that the combined proper motion hc is 
the weighted average of the three proper motions Mi>M2>Mo> 
where we have to use Eq. (26) for the weights of the proper 
motions. 

The mean errors e x .c and e^c of the combined position 
x c{Tc) and combined proper motion hc are derived from Eqs. 
(11), (12), and (29): 



p 2 

£ x,C 



F 2 



1 



Mc Wi + tP2] e x 2 +e 



4 



Mt-T c ) 2 ] c 
1 



-x,2 



1 



«V,i 



''fJ.,2 + Wy.fi 



(31) 



(32) 



-y,l ' c y,2 ~ c fifi 

The weights of xc(Tc) and hc are therefore just the sum of 
the weights of the two positions, 

W X .C = Wx.i + W x .2 , (33) 

and the sum of the weights of the three proper motions, 

Wy.c = Wp,i +«V,2 +^,0 • (34) 
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Due to the separation of the normal equations (21) and (22) 
by our choice of Tc, the values of xc(Tc) an d Mc are not 
correlated. 

In summary, our analytic approach to the combination of 
the two basic catalogues provides a very simple rule: the re- 
sulting combined data are weighted averages of the data in the 
two catalogues, described by Eqs. (19), (23), and (30). For the 
proper motions we have to include the proper motion [i , based 
on the two positions, into this average. The mean errors of the 
combined position and proper motion are given by Eqs. (31) 
and (32). A schematic illustration of the method is shown in 
Fig. 1. While our method is pleasing because of its concep- 
tional simplicity, it is also fully justified by the use of the nor- 
mal equations. As mentioned earlier, the only approximation 
used in our analytic approach is that we have neglected the 
cross-correlations between the coordinate components of the 
positions and proper motions in en* and 5 introduced by HIP- 
PARCOS. 

As a preparation for Sect. 3, we remark (without giving the 
proof here) that our analytic result is identical to the solution of 
another least-square problem in which we treat the catalogue 
data xi(Ti), fi\, x 2 (T 2 ), \i 2 as 'observations'. The new con- 
dition is 

(^i(Ti) - x c {Ti)) 2 + w x , 2 (x 2 (T 2 ) - x c (T 2 )) 2 
+ w n,i (Mi - Vc? + {p-2 - Mc) 2 = min . (35) 

The weights w are given by Eqs. (20) and (26), and xc{Ti) 
follows from 

x c m) = x c (T c ) + ^c(T - T c ) . (36) 

Finally we would like to point out that our analytic method can 
be easily used for a 'decomposition' of a combined catalogue: 
If the catalogues C and 1 are available, then the catalogue 2 can 
be fully reconstructed. In subsequent papers we shall present 
results of this 'inverse combination' method by decomposing 
earlier fundamental catalogues. For example, we have decom- 
posed the FK5 into two parts: the published FK4 and the 'sum- 
mary of subsequent data' which we denote by FK5-4. The FK5- 
4 has also the character of a compilation catalogue. The decom- 
position of the FK5 into a few 'subcatalogues' has the follow- 
ing purpose: If we would like to test the motion of a fundamen- 
tal star for linearity in time or to derive a non-linear motion for 
the star similar to the G solutions of HIPPARCOS, then it is 
helpful to have as many catalogues with different, well-spaced 
epochs as possible at hand. 

3. Numerical approach 

We now present our numerical approach to the problem of com- 
bining two catalogues. In this numerical approach, the correla- 
tions between the five astrometric parameters derived by HIP- 
PARCOS in a standard solution for an apparently single star 
are fully taken into account. 

We change our notation slightly. As a ground-based cata- 
logue, we envisage the FK5, and use the index F (instead of 1). 
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II 
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i 1 1 1 1 1 1 1 

x 2 (T 2 ) 


- 


Ms - 

V( c) ^ 










—f^ Ml 




iii 


i 1 i i i i 1 i i i 



Time t 

Fig. 1. Schematic illustration of the analytic approach to the 
combination of two catalogues. The basic data are the posi- 
tions of a star, x\ and x 2 at the two epochs T\ and T 2 , and 
the corresponding proper motions fi\ and \i 2 . The dotted line 
connects the two positions (open circles). The combined po- 
sition xc (filled circle) at the combined central epoch Tc is 
the weighted average of x\ and x 2 and is located therefore on 
this connecting line. For the positions, error bars are shown. 
The dotted line provides also the direction of the additional 
proper motion ^o> derived from the two positions. The com- 
bined proper motion fic is the weighted average of the three 
proper motions fi\, /i 2 , yu . The dash-dotted line illustrates the 
combined solution xc (t) as a function of time t. 



The HIPPARCOS data receive the index H (instead of 2). The 
index C for the combined catalogue remains. Furthermore we 
will use now the vector and matrix notation in most cases. U is 
the unit matrix with the elements XJu — Ski {k,l — 1, 2, 3, ...); 
is a submatrix filled with zeros. A T denotes the transpose of 
the matrix A, and A -1 its inverse, AA _1 = U. Instead of the 
right ascension a and the corresponding proper motion fi a , we 
use a* = a cos 6 and fi m = [i a cos 8. The parallax of a star 
is denoted by p. 

As already mentioned in Sect. 1, we use in real applications 
only quantities which are differential with respect to HIPPAR- 
COS, e.g. Q!*,F(t) — Q!*,jf(i) instead of a* t p(t). This means, of 
course, that most of the quantities with the index H are zero in 
our scheme. Here we keep, however, the full quantities in order 
not to confuse the reader. 
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The HIPPARCOS Catalogue provides the 5 astrometric pa- 
rameters 



h H (T H ) 



( a* M (T H ) \ 
Sh(T h ) 

Pa*,H 
PH 



(37) 



for the overall epoch T H = 1991.25. Furthermore, the HIP- 
PARCOS Catalogue gives all the data required to derive the 
corresponding variance-covariance matrix D# (Th)- The diag- 
onal line of T>h is occupied by the squares of the mean errors 
of the components of b#(T#), e.g. by 



D H ,i\ — £u*m ( t h) 



(38) 



The non-diagonal elements of D# can be derived from the cor- 
relation coefficients pki,H {k, I = 1, 5) and the mean errors, 
e.g. by 



Dr,12 = Dh,21 — £u*m{Th) £S,h{Th) Pa8M{TH) 

The FK5 provides only 4 astrometric parameters, 

I a. t ^ F (T a . F ) ^ 
Sf{Ts,f) 

Pa*,F 
(J-S,F 



\ 



(39) 



(40) 



since no parallax is given. In contrast to HIPPARCOS, the posi- 
tions are given at the central epochs, T ai p and T$ t F- Hence the 
four FK5 astrometric parameters are not correlated. This means 
that all the non-diagonal elements of the variance-covariance 
matrix Dp are zero. The diagonal elements of T)p are given 
by the squares of the mean errors, e.g. by 



D F ,u = ei^ F (T a . F ) . (41) 
The total vector b of 'observations' (9 elements) is given by 

and the total matrix D (9 x 9 elements) by 
D = 



Dh O 
O Dp 



(43) 



since the HIPPARCOS and FK5 data have no cross-correla- 
tions. The total matrix P of the weights is the inverse of D: 

P = D 1 . (44) 

As the unknowns to be determined, we use the two positions 
at an arbitrary reference epoch, T re f, the two proper motions, 
and the parallax: 

/ a^ C (Tref) \ 
S C'(Tref) 

C = p a *,C • (45) 

M<5,C 

V Pc J 



We now write down the equations of conditions for each 'ob- 
servation', similar to Eq. (4) of our analytic approach. In matrix 
form the equations of conditions are 



b = Ac + v . 

The design matrix A is given by 



(46) 



A = 



(I 







Tref 







o\ 





1 







T H 


— T re f 













1 























1 





















1 


1 





T a ,F 


- T ref 













1 







Tg,F 


— T re f 













1 






















1 


0/ 



(47) 



The vector v contains the residuals of the 9 observational pa- 
rameters. The least-square solution of Eq. (46) requires that the 
sum of the weighted squares of the residuals is minimized: 

v T P v = min . (48) 

From this condition the normal equations follow: 

A T P A c = A T P b . (49) 

The solution c of this equation gives the five astrometric pa- 
rameters of the combined catalogue C for the epoch T re f. The 
mean errors and the correlation coefficients of the elements of 
c are derived from the variance-covariance matrix of c which is 
given by 



Dc-fA^A)- 1 . 
For example, 

£ l\*,c(Tref) = T>CS1 

or 



PaS,C'{T r ef) = D C S2 I (D C ,11 T) C .2l) 



1/2 



(50) 



(51) 



(52) 



The positions, their mean errors and their correlation coeffi- 
cients can, of course, be calculated for other epochs than T re f, 
say for an epoch T, by using 



a*,c{T) = a*. c {T ref ) + p a *,c(T - T ref ) 



(53) 



and a corresponding equation for 5c (T). A specially interest- 
ing epoch is T a ,c,min at which the mean error e a * t c of ct^.c 
reaches its minimum. This epoch corresponds to Tc for a* in 
our analytic approach. In general, T a ^c,min and Ts,c,min are 
not the same epochs. 

The numerical approach is well adapted for many varia- 
tions. Firstly, we can easily add further compilation catalogues 
or individual observational catalogues. Secondly, we can re- 
move unknowns or add further unknowns. For example, the 
'long-term prediction mode' of the FK6 requires that we drop 
the parallax p as an unknown. This can be done either explicitly 
in our equations or by giving pn a weight zero, corresponding 
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FK5 + HIPPARCOS GC + HIPPARCOS 

Quantity a, mean error 5 mean error a, mean error 8 mean error 



Input data: 



Axf(Tf) or Axrr(Trr) 


-7.83 


± 12.51 


+ 96.94 


± 14.78 


-75.60 


± 39.62 


+ 237.68 


± 27.01 


Afi F or Afiac 


+ 0.49 


± 0.40 


- 1.20 


± 0.38 


+ 1.70 


± 0.99 


- 2.83 


± 0.95 


T F or Toe 


1947.84 




1929.73 




1892.60 




1890.30 




AxH{TH,ind) 


0.00 


±0.77 


0.00 


±0.54 


0.00 


± 0.77 


0.00 


±0.54 




0.00 


± 1.01 


0.00 


±0.77 


0.00 


± 1.01 


0.00 


±0.77 




1991.26 




1991.51 




1991.26 




1991.51 






49.48 


±0.99 






49.48 


±0.99 






Results from analytic approach: 


















A^ 


+ 0.18 


±0.29 


- 1.57 


±0.25 


+ 0.77 


±0.40 


-2.35 


±0.27 


Ax c (To) 


-0.03 


±0.77 


+ 0.12 


±0.54 


-0.03 


± 0.77 


+ 0.09 


±0.54 


Afi c 


+ 0.27 


±0.23 


- 1.36 


±0.20 


+ 0.79 


±0.35 


-2.14 


± 0.24 


T C 


1991.10 




1991.44 




1991.22 




1991.47 




Results from numerical approach: 


















Axc(T c ) 


-0.18 


±0.76 


+ 0.14 


±0.54 


-0.27 


± 0.77 


+ 0.15 


±0.54 


Afic 


+ 0.23 


±0.23 


- 1.34 


±0.20 


+ 0.66 


±0.35 


-2.10 


±0.24 


T C 


1991.12 




1991.47 




1991.25 




1991.50 




Ape 


-0.57 


±0.95 






-0.92 


±0.95 







Units: mas, mas/year, or years 



formally to e p .h — * oo. Additional unknowns, e.g. acceleration 
terms similar to the G solutions of HIPPARCOS, are required if 
we would like to determine non-linear solutions for some astro- 
metric binaries. The use of the 'decomposed' FK5 mentioned 
at the end of Sect. 2, demands an increase in both the number 
of catalogues and the number of unknowns. 

4. The effect of cosmic errors 

The combination methods presented in Sects. 2 and 3 are valid 
for single stars only, which move linearily in time on straight 
lines in space. We call this version of our combination method 
the 'single-star mode'. Some binaries can be handled without 
any significant change in this mode of our combination method 
if we use appropriate 'reference points' for the binary, e.g. the 
center-of-mass of the double star. The specific treatment of 
known binaries will be presented elsewhere. 

As mentioned in Sect. 1, undetected binaries introduce cos- 
mic errors into the positions and proper motions, especially into 
the 'instantaneous' HIPPARCOS data. In a subsequent paper 
we shall discuss the modifications of our combination method 
which are required if we take such cosmic errors into account. 
In many applications, it is only necessary to replace the actual 
measuring errors of some quantities by formal errors which in- 
clude both the measuring and the cosmic errors of these quan- 
tities. 



5. An example: a Ari 

In order to illustrate our combination method we give in Ta- 
ble 1 the results for one individual star. All the positions and 
proper motions are given differentially to the HIPPARCOS re- 
sults for the given epoch. Since this is only a numerical conve- 
nience, the mean errors of the quantities given in Table 1 are 
those of the full quantities themselves. The mean errors of the 
ground-based quantities, e.g. of xp, fip, fi , are 'total' values 
which contain both the individual measuring errors and the un- 
certainty in the systematic corrections. The ground-based data 
given are already reduced to the HIPPARCOS system. Hence 
the differences in Table 1 contain only the 'individual' differ- 
ences, not anymore a 'systematic' part. In order to obtain the 
final result of our combination method, one should add the 
results with the index C to the full HIPPARCOS data. e.g. 
xc{Tc) = Ax c (r c ) ± x H (T c ) or [i c = ± \i H . 

The comparison of the results of the analytic approach with 
those of the numerical approach shows a very good agree- 
ment. The numerical approach has made use of the correlation 
coefficients between the HIPPARCOS data for a* (199 1.25), 
5(1991.25), // a «, (J-s, V (+0.20, -0.01, +0.10, +0.26, +0.00, 
-0.35, -0.25, +0.33, +0.02, +0.29). The effect of the corre- 
lations is therefore small. This holds for most stars. The cor- 
relation coefficients for the combined data for a* (199 1.1 2), 
<5(1991.47),/i OT , n s , p (+0.26, 0.00, +0.04, +0.24, +0.04, 0.00, 
-0.16, +0.02, -0.01, +0.08) are smaller than for HIPPARCOS, 
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because the ground-based data are not correlated at all if we use 
the concept of central epochs. 

A comparison of the mean errors of the combined quanti- 
ties with those of HIPPARCOS shows the following: (1) There 
is a significant gain in the accuracy of the proper motion which 
justifies the combination of HIPPARCOS results with ground- 
based data. In our example, the mean errors of the combined 
proper motions are smaller than those of the HIPPARCOS 
proper motion by a factor of 1/4.4 in p a * and of 1/3.9 in ps- 
(2) Since the mean measuring errors of xh are extremely small, 
there is practically no gain in the accuracy of the central posi- 
tions, and we have Tq ~ Th ■ However, for predicting positions 
at epochs which differ from Th by more than a few years, the 
accuracy of the proper motion is governing the error of these 
predicted positions. (3) The improvement of the HIPPARCOS 
parallax pn is here and in most other cases not very signifi- 
cant, since the correlation of pn with the other HIPPARCOS 
data (xh, Ph) is usually small. 

The comparison of the combinations FK5+HIPP and 
GC+HIPP shows the importance of the old observations. The 
FK5 contains essentially all the observations compiled in the 
GC, and in addition many, more accurate, modern observa- 
tions. Nevertheless, the mean errors of the combined data for 
the GC+HIPP are only slightly larger than for FK5+HIPP. The 
analytic approach shows the reason: The accuracy of the proper 
motion p is already very high for the GC+HIPP, mainly be- 
cause of the large epoch difference of the central positions. 



nation of the FK5 data with the HIPPARCOS observations, are 
typically by a factor of about two more accurate than the proper 
motions in either the FK5 alone or in HIPPARCOS Catalogue 
alone. Already p is typically slightly more accurate than pn- 



Table 2. Error budget for FK6 proper motions in the 'single- 
star mode' 

Typical mean errors of proper motions 
(in one component, averaged over p at and pg', units: mas/year) 



Sample of stars: 



1535 FK 



1202 FK 



rms aver. median rms aver. median 



HIPPARCOS 



FK5 



Mo 



random 
system 
total 

random 
system 
total 



FK6 = FK5+HIP 



0.82 

0.76 
0.28 
0.81 

0.53 
0.24 
0.58 

0.35 



0.63 

0.64 
0.25 
0.70 

0.43 
0.23 
0.49 

0.33 



0.68 

0.77 
0.28 
0.83 

0.54 
0.25 
0.59 

0.35 



0.61 

0.67 
0.25 
0.72 

0.45 
0.23 
0.51 

0.34 



6. Error budget 

In order to show the overall improvement in the accuracy of the 
proper motions obtained by using our combination method, we 
present here the statistics of the appropriate mean errors for a 
few samples of stars. 

The mean errors e M of the proper motions p, given in Tables 
2 and 3, refer to one 'mean' coordinate component of p. For 
each star, the 'mean' value e^m is obtained as a root-mean- 
square (rms) average over e^a* and e^^. Then these individ- 
ual values of e^m are averaged over the sample of stars un- 
der consideration, either by taking an rms average over all the 
£^,id or by selecting the median value of e^m in the sample. 

In Table 2, we present the error budget for the combination 
of the FK5 with HIPPARCOS for the basic FK5 stars. This 
corresponds to the results of the 'single-star mode' of the FK6 
(Wielen et al. 1998). We consider two samples of stars: (a) all 
the 1535 basic FK5 stars, and (b) 1202 'apparently single' basic 
FK5 stars. The error budget of all the basic FK5 stars is slightly 
fictious in so far as this sample also contains binaries which 
cannot be safely treated by a direct combination method. Nev- 
ertheless, the error budget for the 1535 stars is still a valid indi- 
cator for the overall accuracy of our combination method. The 
error budget for the 1202 stars is slightly biased in favour of the 
HIPPARCOS Catalogue, since most of the stars for which es- 
pecially HIPPARCOS is less accurate than usually (e.g., C so- 
lutions for visual binaries) are removed from this sample. Table 
2 shows that the FK6 proper motions, obtained from the combi- 



ratio of HIPPARCOS 2.3 1.9 1.9 1.8 

to FK6 errors 



In Table 3, the error budget for the combination GC+HIP 
of the GC data with the HIPPARCOS observations are shown. 
The samples of stars are the following: (a) all the 29 717 GC 
stars which we have identified in the HIPPARCOS Catalogue; 
(b) 1 1 773 stars in the GC which have standard solutions by 
HIPPARCOS (i.e., they do not occur in the DMS Annex of 
the HIPPARCOS Catalogue, which would indicate a double or 
multiple object) and which are not known to be double from 
ground-based observations and which have accurate values of 
Po(GC) (mean errors smaller than 2.0 mas/year in each coor- 
dinate); (c) all the 1534 basic FK5 stars in the GC (one ba- 
sic FK5 star is not in the GC); (d) 1201 'apparently single' 
basic FK5 stars in the GC. Table 3 shows that for the 1 1 773 
apparently single GC stars the combination GC+HIP produces 
proper motions which are more accurate than the HIPPARCOS 
proper motions by about 20 %. The overwhelming part of this 
improvement stems from the proper motion Pq(gc)- F° r me 
brighter GC stars, represented by the FK stars in Table 3, the 
gain in accuracy is much larger, typically by a factor 1 .5 to 2 . 
A comparison of the results for GC+HIP given in Table 3 with 
those for FK6 = FK5+HIP in Table 2 for the basic FK5 stars 
confirms that the older ground-based observations carry a very 
large weight in the combination with the HIPPARCOS data: 
For the basic FK5 stars, the mean errors of the proper motions 
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Typical mean errors of proper motions 
(in one component, averaged over and \ units: mas/year) 



Sample of stars: 
Median of m v : 



HIPPARCOS 
GC 

random 
system 
total 



29717 GC 
7 m 



11773 GC 

6 m 8 



1534 FK from GC 
4™8 



1201 FK from GC 

4« l 9 



rms average median rms average median rms average median rms average median 



1.47 

10.57 
0.38 
10.58 



0.73 

9.38 
0.34 
9.39 



0.75 

8.59 
0.38 
8.60 



0.69 

7.55 
0.33 
7.56 



0.82 

3.33 
0.47 
3.37 



0.63 

2.31 
0.40 
2.38 



0.68 

3.44 
0.47 
3.47 



0.61 

2.45 
0.40 
2.48 



random 
system 
total 



1.78 
0.15 
1.78 



1.75 
0.11 
1.75 



1.43 
0.12 
1.43 



1.45 
0.11 
1.45 



0.66 
0.19 
0.66 



0.55 
0.14 
0.55 



0.67 
0.18 
0.67 



0.57 
0.14 
0.57 



GC+HIP 



0.72 



0.63 



0.62 



0.58 



0.42 



0.39 



0.41 



0.39 



ratio of HIPPARCOS 
to GC+HIP errors 



2.0 



1.16 



1.21 



1.19 



2.0 



1.6 



1.7 



1.6 



of the combination GC+HIP are only slightly larger (by about 
20 %) than those of the FK6 = FK5+HIP. 

The error budgets for the stars in the bright and faint ex- 
tension of the FK5 and for the remaining FK4Sup stars are not 
given here. For these stars the gain in accuracy of their proper 
motions by using our combination method lies inbetween the 
gain for the basic FK5 stars and that for the GC stars. 

We would like to warn those readers who would like to test 
our combination method with the help of the numbers given 
in Tables 2 and 3, using the equations given in Sect. 2 . The 
quantities given in Tables 2 and 3 are ensemble averages of 
the mean errors e M , while our combination method works indi- 
vidually star by star. Since < > differs in general from 
1/ < >, such a test is not meaningful. A numerical check of 
our method can, of course, be carried out by using the numbers 
given in Table 1 for the individual star a An. 

It should be emphazised that the gain in accuracy is sig- 
nificantly larger in the 'long-term prediction mode' (Wielen et 
al. 1998) in which we take into account the 'cosmic errors' in 
the HIPPARCOS results (see Sect. 4). The error budget for the 
long-term predictions will be discussed in detail in a subse- 
quent paper. 

7. Conclusions and outlook 

We have shown that the combination of ground-based compila- 
tion catalogues with the HIPPARCOS data produces for many 
stars proper motions which are significantly more accurate than 
the HIPPARCOS proper motions themselves. For stars which 
are bright (~ 5 m ), well-observed and truly single, the gain in 



accuracy is about a factor of two. For stars of intermediate mag- 
nitudes (<~ 7 m ), the smaller gain of about 20 % on average is 
still valuable for many applications. 

The combination method described above has been used 
intensively and successfully in the construction of the FK6, 
in which the basic FK5 is combined with HIPPARCOS. The 
method described here provides the FK6 results in the 'single- 
star mode' (see Wielen et al. 1998). The first part of the FK6 
will be published in the near future. The second part of the FK6 
deals with double stars and requires a slightly modified version 
of our combination method. 

In order to take into account the cosmic errors in the HIP- 
PARCOS results, which are caused by undetected astrometric 
binaries, our combination method has to be changed slightly. 
The corresponding procedures, which lead to the 'long-term' 
and 'short-term' prediction modes (Wielen et al. 1998, 1999), 
will be presented in a subsequent paper. 

References 

Bien R., Fricke W., Lederle T., Schwan H., 1978, Veroff. Astron. 

Rechen-Inst. Heidelberg No. 29 
Boss B., Albrecht S., Jenkins H., Raymond H., Roy A.J., Varnum 

W.B., Wilson R.E., 1937, General Catalogue of 33 342 Stars for 

the Epoch 1950, Carnegie Institution of Washington, Publ. No. 

486 

Eichhorn H., 1974, Astronomy of Star Positions, Frederick Ungar 

Publ. Co., New York 
ESA 1997, The Hipparcos Catalogue, ESA SP-1200 
Fricke W., Schwan H., Lederle T., Bastian U., Bien R., Burkhardt 

G, du Mont B., Hering R., Jahrling R., JahreiB H., Roser S., 



10 



R. Wielen et al.: Combination of ground-based astrometric compilation catalogues with the HIPPARCOS Catalogue 



Schwerdtfeger H.M., Walter H.G., 1988, Veroff. Astron. Rechen- 
Inst. Heidelberg No. 32 
Fricke W., Schwan H., Corbin T., Bastian U., Bien R., Cole C, Jack- 
son E., Jahrling R., Jahreifi H, Lederle T., Roser S., 1991, Veroff. 
Astron. Rechen-Inst. Heidelberg No. 33 
Kopff A., 1937, Veroff. Astron. Rechen-Inst. Berlin No. 54 
Kopff A., 1938, Abh. PreuB. Akad. Wiss., Jahrgang 1938, Phys.-math. 
Klasse, No. 3 

Kopff A., Nowacki H., Strobel W., 1964, Veroff. Astron. Rechen-Inst. 

Heidelberg No. 14 
Schwan H, Bastian U., Bien R., Jahrling R., Jahreifi H., Roser S., 

1993, Veroff. Astron. Rechen-Inst. Heidelberg No. 34 
Wielen R., 1988, In: IAU Symposium No. 133, Mapping the Sky, 

eds. Debarbat S., Eddy J.A., Eichhorn H.K., Upgren A.R., Kluwer 

Publ. Comp., Dordrecht, p. 239 
Wielen R., 1995a, A&A 302, 613 

Wielen R., 1995b, In: Future Possibilities for Astrometry in Space, 

eds. Perryman M.A.C., van Leeuwen F, ESA SP-379, p. 65 
Wielen R., 1997, A&A 325, 367 

Wielen R., 1998, Tatigkeitsbericht des Astron. Rechen-Inst. Heidel- 
berg fur das Jahr 1997 = Mitt. Astron. Ges. No. 81, 327 

Wielen R„ Schwan H., Dettbarn C, Jahreifi H., Lenhardt H., 1997, In: 
Hipparcos Venice '97, Presentation of the Hipparcos and Tycho 
Catalogues and first astrophysical results of the Hipparcos space 
astrometry mission, eds. Battrick B., Perryman M.A.C., Bernacca 
PL., ESA SP-402, p. 727 

Wielen R., Schwan H., Dettbarn C, Jahreifi H., Lenhardt H, 1998, In: 
The Message of the Angles - Astrometry from 1798 to 1998, Pro- 
ceedings of the International Spring Meeting of the Astronomi- 
sche Gesellschaft, held in Gotha, 11-15 May 1998, eds. Brosche 
P., Dick W.R., Schwarz O., Wielen R., Acta Historica Astrono- 
miae, Vol. 3, Verlag Harri Deutsch, Thun and Frankfurt am Main, 
p. 123 

Wielen R., Schwan H., Dettbarn C, Jahreifi H., Lenhardt H, 1999, 
In: Modern Astrometry and Astrodynamics, Proceedings of the 
International Conference honouring Heinrich Eichhorn, held at 
Vienna Observatory, Austria, 25-26 May 1998, eds. Dvorak R., 
Haupt H.F, Wodnar K., Verlag der Osterreichischen Akademie 
der Wissenschaften, Wien, p. 161 



